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We present experimental results on the preparation of a desired quantum state in a two-level 
system with the maximum possible fidelity using driving protocols ranging from generalizations of 
the linear Landau-Zener protocol to transitionless driving protocols that ensure perfect following of 
the instantaneous adiabatic ground state. We also study the minimum time needed to achieve a 
target fidelity and explore and compare the robustness of some of the protocols against parameter 
variations simulating a possible experimental uncertainty. In our experiments, we realize a two- 
level model system using Bose-Einstein condensates inside optical lattices, but the results of our 
investigation should hold for any quantum system that can be approximated by a two-level system. 
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I. INTRODUCTION 

The "simplest non-simple quantum problem" is the 
evolution of a two-level quantum system with a time- 
dependent Hamiltonian [1 . A paradigmatic example of 
such a system is the Landau-Zener problem, in which two 
levels whose energy depends on time through some pa- 
rameter experience an avoided crossing during the tem- 
poral evolution. If the time variation of the energy levels 
is linear, the resulting evolution gives rise to the well- 
known phenomenon of Landau-Zener (LZ) tunneling |2J- 
4, a scenario lying at the heart of a host of quantum 
phenomena. The LZ tunneling probability is often used 
to estimate the degree of adiabaticity achievable when 
a two-level system is forced through an avoided crossing 
of its energy levels. In modern terminology, this can be 
described as a problem in quantum control: if the goal 
is to keep the system as close as possible to the low- 
est energy level before and after the crossing, one has to 
choose a time dependence that minimizes the probability 
of LZ tunneling. In other words, one aims to maximize 
the fidelity of the control protocol, defined as the over- 
lap between the desired and the actual final state. For 
the standard LZ problem with a linear time-dependence 
of the energy, the goal of perfect preparation of the fi- 
nal state can only be achieved for infinitely slow sweep 
speeds or for infinitely large coupling strengths, both of 
which are impractical (and, indeed, unphysical). 

This observation leads naturally to the following ques- 
tion: can the degree of fidelity in a LZ type quantum 
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control protocol be improved by using non-linear sweeps 
in some control parameter T? If the off-diagonal element 
lo of the Hamiltonian can be also varied during the proto- 
col, what is the optimum combination of the Hamiltonian 
parameters that maximizes the fidelity of the final state? 
Finally, is it possible to choose the parameters of the two- 
level system in order to speed up the transfer from the 
initial state to the final one? All these questions repre- 
sent key issues within the present effort to control quan- 
tum phenomena in different areas from atomic physics to 
solid-state physics and molecular physics [SJ [7] . 

The quest for optimal control of the two-level systems 
has a long history, and protocols denoted as Rosen-Zener, 
Allcn-Ebcrly or Demkov-Kunike have been studied for 
many years |S]. In magnetic resonance, composite pulses 
were introduced to create fast and robust transfer pro- 
tocols Within the framework of adiabatic quantum 
computation, Roland and Cerf |10) introduced an adia- 
batic protocol that is much faster than the standard LZ 
protocol. Protocols with nonlinear time dependence of 
the r parameter, the so-called generalized LZ sweeps, 
have been investigated because, in contrast to the LZ 
problem with a linear sweep, they alllow a complete tran- 
sition to the target state with a finite sweep rate [TT| [T2] . 
An alternative approach is to use shortcut-to-adiabaticity 
protocols leading to the same target state as the LZ 
approach, but in a shorter time. Several methods to 
find shortcuts to adiabaticity have been proposed for 
two-level systems, such as reverse engineering using the 
Lewis-Riesenfeld invariants [T71 [T5] . Finally, the super- 
adiabatic (also known as counter-adiabatic or transition- 
less) protocols for a given time- varying Hamiltonian con- 
struct an auxiliary Hamiltonian that exactly cancels the 
non-adiabatic part of the original Hamiltonian and thus 
ensures perfect adiabatic following [13HT6] . 

A related issue is the desire to (accurately) realize these 
quantum manipulations in the shortest possible time. 
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There is, however, a fundamental limit towards this goal 
which is set by quantum mechanics. There is a maximum 
speed, imposed by the Schrodinger evolution, at which a 
quantum system can evolve |19H23j . If a quantum sys- 
tem can be controlled by a set of external parameters, 
the speed limit can be attained only for a suitably tai- 
lored time-dependence of these control parameters. 

All the above questions have so far received limited 
experimental attention. The Roland-Cerf protocol was 
experimentally implemented on different quantum com- 
putation algorithms by nuclear magnetic resonance tech- 
niques in |24j . The control of nonadiabatic transitions in 
an adiabatic transformation was performed in the exper- 
iment of ref. |25j . A nonadiabatic transformation lead- 
ing to a nearly perfect preparation of the target state 
was realized in an experiment on a Bose-Einstein con- 
densate within an optical lattice in [26] . Our recent work 
with Bosc-Einstein condensates compared the efficiency 
reached in the LZ, Roland-Cerf and superadiabatic pro- 
tocols 27J. The shortcuts to adiabaticity based on the 
the Lewis-Riesenfeld invariants for the case of Hamilto- 
nian with a dependence also on the spatial coordinates 
were investigated in experiments on trapped ultracold 
gases [2"5] . 

The quantum speed limit had not been explicitly tested 
before our investigation in ref. [2 7) . In nuclear magnetic 
resonance or atomic physics, it is typically assumed that 
the most rapid transfer is produced by a n Rabi pulse. 
In fact, for a two- level Hamiltonian with constant pa- 
rameters, the 7r pulse corresponds to the quantum speed 
limit. However, for more complicated protocols such as 
those investigated here, the deviation from the quantum 
speed limit is an important issue. 

The present work reports an experimental investiga- 
tion of different protocols regarding their fidelity and 
the deviation from the quantum speed limit. Section II 
briefly introduces the quantum driving of our two level 
paradigmatic quantum system and the also the quantum 
speed limit. Starting from the linear LZ as a reference 
protocol investigated in [29j [30] , we explored generalized 
LZ protocols where lo is fixed and T is scanned non lin- 
early (Section III). Section IV discusses the Roland and 
Cerf protocol. Section V introduces the idea of supera- 
diabatic protocols and reports time-dependent measure- 
ments for both the diabatic and adiabatic bases. Section 
VI investigates non-optimized Roland and Cerf protocols 
and superadiabatic protocols, for a deviation of control 
parameters from optimal values. Section VII discusses 
the best use of resources, which are always limited ex- 
perimentally, in terms of the minimum required to reach 
a certain target fidelity. Section VIII presents the conclu- 
sions of our experimental and theoretical investigations. 



II. THE SYSTEM AND THE TARGET 

A. Hamiltonian 

As described in the Appendix and shown schematically 
in Fig. [T] our system is realized with ultracold atoms 
forming a Bose Einstein condensate (BEC), in an accel- 
erated optical lattice. Under appropriate conditions the 
BEC can be described as a two-level system defined by 
the two lowest Bloch bands: 

U = hT(t)a z + huj{t)a x , (1) 

where F is the energy between the diabatic components 
of the Hamiltonian and u) is the intensity of the cou- 
pling between the diabatic components, causing the 
splitting between the bands. The above generalized 
time-dependent LZ Hamiltonian has instantaneous 
eigenvectors \ip giC (t)) , its adiabatic eigenstates. The 
diabatic states |0), |1) are the eigenvectors of the Hamil- 
tonian with uj(t) = 0. 

The case of constant lo and T(t) depending linearly on 
t corresponds to the standard LZ Hamiltonian [H [21 [5] 
with an avoided crossing with a minimum energy gap of 
2fko at T(t) = 0. For our system of a BEC in an optical 
lattice, the lattice depth controls the tunneling barrier 
and hence the gap, while an acceleration of the lattice 
controls the time dependence of T. 

Having prepared the system in the ground state of % 
for a given choice of the initial parameters and we 
would like to reach, after an evolution of duration T, 
the ground state of H for the final parameters LOf and 
Tf. The initial and final values of the parameters are 
chosen to be on opposite sides of the anticrossing. The 
different protocols examined in the following correspond 
to different time dependencies T(t) and ui(t). 

For simplicity, T, lo and T are dimensionless, using 
the natural units of energy and time in our Rubidium 
system, see Appendix A. In addition, for a given value 
of T, we introduce the rescaled time r = t/T, r € [0, 1], 
and we assume that the sweep ranges from —Tq to Fo, 
with r = 2, and is point-symmetric about t = 0.5. 



B. Fidelity 

We characterize the protocol efficiency for the system 
going from the initial state l^m^r = 0)) = \ip g (r = 0)) 
to the final state | x I'fi n (r = I)) through the following final 
fidelity: 

^nH(^(T=l)|* fi „(T=I))| 2 . (2) 

The evolution of the wavefunction along the adiabatic 
trajectory at an intermediate time r for given T may be 
tested by extending the above definition to a time depen- 
dent fidelity Fir) . In the experiment the wavefunction at 
time r is determined by a projective measurement on the 
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FIG. 1: Experimental realization of an effective two level sys- 
tem exploiting the band structure of a Bose-Einstein conden- 
sate in an optical lattice. The avoided crossing is located at 
the edge of the Brillouin zone and the two adiabatic states 
|i/)g ie (r)} are the lowest energy bands. The initial and final 
states of the time evolution | ^/ini.nn) are also indicated. 
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FIG. 2: Time dependence of the control parameters for differ- 
ent protocols. In (a) F(r) for power law sweeps: red dotted 
line a — 2, blue continuous line a — 3, green dot-dashed line 
a — 4. In (b) F(r) for linear+sin sweeps: red dotted line 
5 — 0.1, blue continuous line S = 0.2, green dot-dashed line 
5 — 0.4. The black dashed line corresponds to the linear LZ 
sweep. 



adiabatic or diabatic bases of the Hamiltonian of Eq. ([I]). 
Notice that for the case of the diabatic basis, the projec- 
tive measurement of the wavefuction \ip g ) results in the 
following probability: 



->diab 



(r) = |(l|^ g (r))| i 



(3) 



where we denote by |0) the state that at t = has 
the lowest energy. For the superadiatic protocols the 
wavefunction follows precisely the instantaneous adia- 
batic eigenvalues of the Hamiltonian of Eq. (JlJ, so that 
the above probability becomes 



■P diab (r) = 



r 2 (r) + w 2 (r) + r(r)Vr 2 (r)+w 2 (r) 



(4) 



by two short pulses (in theory delta-functions) with a 
pulse area of 7r/4. The transfer time T associated with 
the l^ini) to | 'I' fi n ) transfer through the above protocol is 



T = 2t, 



o 



arccos|(W nn |*i 



(6) 



A Rabi rotation of an angle tt corresponds to a time 
Tn = tt/lo, and in fact in the limit of to — > 0, the above 
time coincides with T v . For uj values up to 1, the T v 
dependence on u> matches the lower continuous line of 
Fig. 10 At larger w values, \tp g (T = 0)) differs from |0), 
and the initial tilting angle 9 = lu/Tq of the Bloch vector 
leads to a required Rabi rotation smaller than 7r. 



C. Quantum speed limit 



III. GENERALIZED LZ SWEEPS 



The quantum speed limit (QSL) protocol allows us to 
take the system from the initial state to the final state, 
with final fidelity J-" nn = 1, in the shortest possible time 
T. This lower bound is rooted in the Heisenberg uncer- 
tainty principle [TM23"] . For the case of constant to, we 
found that the protocol minimizing T is [27 



r( T ) = 




for t = 
for t G [0,t ] 
forte [t ,T-t ] 
forte [T-t ,l] 
for t = T , 



(5) 



We begin our investigation of the dynamics of generic 
two-level systems starting from the simplest possible gen- 
eralizations of the LZ scenario. In this Section we keep 
the coupling ui constant and study deviations of T(t) from 
the linear dependence assumed in the original LZ prob- 
lem. We introduce nonlinear sweeps in two ways: by 
making the functional form of T(r) follow some power 
law in r and by adding a small nonlinear contribution to 
the original linear sweep. 



Power laws 



where Tm and t are, respectively, asymptotically large 
and small quantities which satisfy the condition r^to = 
tt/4. 

This 'composite pulse' protocol, in close analogy to 
composite pulses in NMR [5] , represents half a Rabi oscil- 
lation with frequency u> and T=0, preceded and followed 



The first scheme leads to sweeps of the following form 
satisfying the boundary conditions T(0) = —2 and T(l) = 
2: 



r(r) = 



-2 q+1 (1/2-t) q for0<r<l/2 
2 q+1 (t-1/2) q forl/2<T<l, 



(7) 
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The case a = 1 corresponds to the linear LZ sweep and 
different power laws lead to the time dependencies de- 
picted in Fig. [2ja) . For a > 1 the bandgap at the edge 
of the Brillouin zone is crossed once with zero speed, and 
for increasing values of a the width of the low speed re- 
gion increases, i.e. |r(r)| < C, for \t - 1/2| < (C/2) 1 / a . 
For example, if a = 4 and C = 10~ 3 , the width of the 
low speed region is 15% of the total sweep duration. 

Power law sweeps have been studied theoretically by 
Garanin and Shilling [IIJ [12] . They calculated the cor- 
rections to the LZ transition probability due to the non- 
linearity, finding that it has an oscillatory behavior if the 
sweeps exhibits a significant retardation in the vicinity of 
the resonance, i.e. T(t) is close to zero around t = 0.5, 
in agreement with our experimental observations. 

We have examined experimentally and numerically the 
dynamics of a two level system for various values of a. 
In Fig. [3ja) the fidelity of the final state for a = 1, 2, 4 
is plotted as a function of the total duration T of the 
sweep, keeping u) — 0.5 fixed. Numerical simulations of 
the power laws sweeps shown in Fig.|3ja) agree well with 
the experimental data. Sweeps with larger values of a 
reach a given final fidelity in a shorter time. Fixing a 
threshold value at J-fi n = 0.9, the required sweep dura- 
tion T . 9 strongly depends on a, as shown in Fig. [3jb) . 
For large values of a, T g approaches the minimum time 
given by the quantum speed limit for the parameters of 
our system (Tqsl — 2.75). In that regime, the protocol 
approaches the form of the composite pulse protocol of 
Eq.0 



B. Linear+sin 



Nonlinear sweeps can also be created by adding a small 
sinusoidal contribution to the original linear sweep: 
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FIG. 3: (a) Final fidelity as a function of the duration T of 
the power law sweep. Blue triangles a = 4, red circles a — 2, 
black squares for linear LZ with a = 1. The continuous lines 
(blue, red and black) represent numerical simulations for the 
corresponding power law sweeps, (b) Minimum time To. 9 to 
achieve J-fi n = 0.9 for the power law sweeps. The experimental 
points (solid circles) are compared to a numerical simulation 
(dot-dashed red line). The black line shows the minimum 
time necessary to achieve J-fin = 1 according to the quantum 
speed limit. All data with cj = 0.5 



r(r) = 4[r + (5sin(27rr)] - 2 . 



(8) 



IV. ROLAND-CERF PROTOCOL 



which satisfies the requirements T(0) = — r(l) = 2 and 
r(0.5) = 0, as shown in Fig. [2jb) . Depending on the 
value of 8 the eigenenergies of the Hamiltonian of Eq. ([I]) 
are characterized by either one or three avoided crossings 
during the temporal evolution. In our experiment the 
bandgap is crossed only once for S = 0.1 and twice for 
6 = 0.2 and 0.4. The speed of the crossing at r = 0.5 
depends on S as T(0.5) = 4 — 2n5. Our measurements 
and the numerical simulation of Fig. show that the 
fidelity of the final state is not monotonic as a function of 
the duration T of the sweep. In addition, for the largest 
S the fidelity experiences an oscillatory behavior similar 
to that of power law sweeps. Those oscillations allow us 
to reach the final state more rapidly than the linear LZ 
sweep. 



The locally adiabatic protocol proposed by Roland and 
Cerf [10] in the context of adiabatic quantum computa- 
tion poses stricter constraints on the time evolution of 
the system. We summarize the analysis in |10| adapting 
it to the case of the Hamiltonian H of Eq. [l] under the 
assumption that the coupling oj is kept constant. At any 
instant of the evolution the fidelity of the state with the 
instantaneous ground state is required to be 

Hr)^\mr)\Ur))\ 2 = l~e 2 , (9) 

It follows [57] that the total time of the protocol is 
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the realization of the Roland and Cerf (RC) protocol 
requires a precise knowledge of the parameters of the 
underlying physical system, since the duration of the 
sweep Tp has to be matched to the coupling constant u. 
Moreover, we found that, for a given uj, the RC protocol 
reaches the target fidelity T = 0.9 within a time only 
twice that of the time-optimal composite pulse protocol. 
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FIG. 4: In (a) numerical simulations for the final fidelity 
as a function of the duration T of the linear plus sin sweep, 
for different values of the parameter 8 (8 = 0.1 red dash- 
dotted line, 8 = 0.2 blue dashed line, 8 = 0.4 green continuous 
line). For increasing values of 8 the final state is reached at 
an earlier time. In (b-d) experimental results are compared to 
the simulations. In (b) 8 = 0.1, in (c) 8 = 0.2, in (d) 8 = 0.4. 
Typical error bars are reported. The dashed black line shows 
the theoretical result for the linear LZ sweep. All data with 
lj = 0.5. 




FIG. 5: Time dependence of the control parameters for 
Roland-Cerf protocol: blue rj 2 — (4 + a; 2 ) -1 corresponding 
to the optimized one, red rf = 0.85(4 + u 2 ) -1 , dashed-dotted 
green rj 2 = 1.05825(4 + a; 2 )" 1 » VO^E- The black dashed 
line corresponds to the linear LZ sweep. 



for a time-dependence of T(t) of the form 



T(r) 



Wl>(r - 1/2) 



(11) 



y/1 - 16e 2 w 2 T£(T- 1/2) 2 ' 
We have implemented the above protocol and verified 



V. SUPERADIABATIC PROTOCOLS 

A further generalization of the LZ scenario is realized 
by allowing the possibility of a time- varying Li(t). We 
have analysed those protocols that ensure a perfect fol- 
lowing of the adiabatic ground state |V'g(' r )) f° r au r i 
called 'superadiabatic' or transitionless with reference to 
Berry's work jT3J Q3] and counter-adiadiabatic with ref- 
erence to the work of Demirplak and Rice [IS] US] ■ Such 
protocols enable fully adiabatic preparation of a desired 
quantum state of a system in a finite time. The under- 
lying idea is that, given a Hamiltonian %, it is always 
possible to construct an additional control Hamiltonian 
V. c such that the under the action of % + H c the sys- 
tem will remain perfectly in the ground state of H at all 
times. In principle, H c should be experimentally imple- 
mented through an extra control field. In our previous 
work [27], we found that it is possible to eliminate the 
need for an additional field by suitably transforming the 
time dependence of the initial Hamiltonian r — > V and 
u) — > ui'. 

Two superadiabatic protocols are considered here. 
The superadiabatic linear protocol, corresponding to the 
transformation of the linear LZ protocol with T(r) = 
4 (r — 1/2), is given by 

r'(r) =r(r)- T2[(T T )2+ \^ +1 , 

<A"r) =^1+t^ t -iV+^) 2 - (12) 

The superadiabatic tangent protocol, corresponding to 
the transformation of the tangent protocol for which T(t) 
is not modified by the transformation, is given by 

1 



r'(r) = r(r) = wtan 2 t arctan 



arctan(-) 2 

Ma) ' 

{Tuf 



(13) 



For both of these protocols the discontinuities in the 
derivative of T at r=0 and r=l lead to delta- functions in 
r' at the beginning and at the end of the time evolution, 
which can be realized in practice using large but finite 
corrections ATm for a short duration At such that 



the Eq. (10 1 for the Tj- dependence on w [57] ■ However, 



At Ar m = T 2 arc tan 



Tuj — ujT 
2w(r 2 + u 2 



(14) 
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FIG. 6: Time dependence of r'(r) (black) and u/(r) (gray) 
for the superadiabatic linear (dashed line) and the superadi- 
abatic tangent protocols (continuous line). 



where the - and + signs refer to the correction at the 
beginning and at the end of the protocol, respectively. 
These superadiabatic protocols are shown schematically 
in Fig. [I 



A. Diabatic and adiabatic bases 

The dynamics of the system defined by the Hamilto- 
man of Eq. can be measured in the diabatic and adi- 
abatic bases. The dynamical evolution may be frozen at 
different times by performing a projective quantum mea- 
surement on the states of a given basis, as we explored in 
ref. [H] for our implementation of a two-level quantum 
system based on BECs in an optical lattice. We have 
measured the fidelity JF{t) during the time evolution of 
the protocol in both bases. 

For measurements in the adiabatic basis, after BEC 
preparation and acceleration, the lattice depth was adi- 
abatically reduced to zero, within 400 (is (7.9 in natural 
units). The protocol fidelity J-(t) was obtained as the 
fraction of atoms remaining in the lowest Bloch band, 
which corresponds to the fraction of atoms in the q = 
momentum class relative to the total atom number. Re- 
sults for the linear LZ protocol and the superadiabatic 
protocols are shown in in Fig. [7ja). 

For the measurement in the diabatic basis, after the 
initial loading phase and acceleration up to the time r, 
the atomic sample was projected onto the free-particle 
diabatic basis by switching off the optical lattice instan- 
taneously (within less than 1 fis, corresponding to 0.02 
in natural units). The number of atoms in the q = 
and q = 2p rec momentum classes were measured, and 
the probability 7? dlab (r) derived as the fraction of atoms 
in the q = 2p rec velocity class relative to the total atom 
number. Diabatic results for the linear LZ protocol and 
the two superadiabatic ones are shown in Fig.^b). The 
measured values are in good agreement with the predic- 



c 0.8 - 



0.6 - 




FIG. 7: (a) Fidelity vs r during the protocol for the superadi- 
abatic linear (red circles), superadiabatic tangent (green tri- 
angles) and linear LZ (black squares) protocols in the adia- 
batic base. p dlab as a function of the time, in (b), and of 
the condensate quasimomentum q, in (c), for the superadi- 
abatic linear (red circles) and superadiabatic tangent (green 
triangles) protocols, the continuous lines derived from Eq. Q . 
The black squares and lines show the reference LZ protocol. 
Error bars are one standard deviation. Data for cj = 0.5 and 
T = 5.9. 



tions based on Eq. Q. 

In Fig-J^c) we show the experimental data for •p dlab (-r) 
during the evolution of q at time r, and compare them 
to a theoretical prediction based on Eq. Q and on the 



linear relation between q(r) and T(r) given by Eq. (A2|. 
Since to is constant for the LZ and the tangent proto- 
cols, the states evolving under different superadiabatic 
Hamiltonians have the same parametrization in q. 



7 




FIG. 8: (Color online) (a) Measured final fidelity vs T for 
the Roland-Cerf protocol with uj — 0.5 and different values 
of 77. Green triangles if = 0.1, black diamonds 77 2 — 0.2, red 
circles 77 = ?7opt, blue squares rj 2 = 0.249. The continuous lines 
report the theoretical predictions, with the best fit obtained 
supposing ui = 0.45. The T — 0.9 target, denoted by the 
dashed line, is reached at shorter times T increasing 77, with 
an oscillating J-n n (T) occurring for 77 > rjopt = V0. 237954. 
(b) To. 9 vs rf . The limiting value of the optimized protocol 
for the given ui is denoted by dashed red line, with the Tsql 
limit denoted by the dash-dotted black line. 



VI. NON-OPTIMIZED PROTOCOLS 



As in realistic experimental conditions parameters may 
be known only to within some uncertainty, protocol ro- 
bustness with respect to parameter variations is relevant. 
The RC protocol and the superadiabatic protocols re- 
quire the knowledge of the coupling constant uj in or- 
der to compute the parameter sweeps. We can simulate 
a possible experimental uncertainty by deliberately per- 
forming the protocol with parameters deviating from the 
optimal values. 



A. Roland-Cerf protocol 

The experimental implementation of the RC protocol 
demonstrated the high sensitivity of the final result on 
the protocol parameters. In order to characterize our 
analysis, we introduce the following parameter appearing 
in eq. [TT] 



77 = euTjr, 

leading to the following time dependence 



r(r) = 



4y/l -4?7 2 (t- 1/2) 
y/l - 16t7 2 (t- 1/2) 2 



(15) 



(16) 



Notice that the above T(t) dependence satisfies the 
boundary conditions T(0) = — r(l) = 2 for any value of 
77. The value rj op t corresponding to the optimized Roland- 
Cerf protocol can be expressed as a function of u> only: 



1 



r\ ov t = 



(17) 



For 77 ^= r/opt, Eq- (16) defines a non-optimized Roland- 
Cerf protocol, with 77^ < 0.25, as it may happen if uj is 
not precisely known. 

Fig.[2]Jc) reports the T(r) dependence for different val- 
ues of 77, including the one corresponding to the optimized 
protocol. Notice the close resemblance between the F(t) 
sweep for 77 V0.25 and the composite pulse protocol of 
Eq. ([5] 

Fig. Ma) shows the experimental results of the final fi- 
delity .Ffi n as a function of the duration T of the protocol 
for different values of 77. It is clear that for increasing 77 
the J-fi n = 0.9 threshold is reached at decreasing values 
of T. At 77 larger than rj opt the threshold is reached for 
smaller T, but in that regime there are strong oscillations 
in the final fidelity. The measured dependence of To. 9 on 
?7 2 is reported in Fig[8^b) and compared to the results of 
a numerical simulation. 



B. Superadiabatic protocol 

For comparison with the Roland-Cerf case we have also 
tested the system time evolution and Ta n for the su- 
peradiabatic protocols for the case where the parameters 
of the superadiabatic Hamiltonian deviate from those of 
Eqs. (12) and (13). In order to test the sensitivity of 



the superadiabatic tangent protocol to a variation in the 
control parameters, we measured the final fidelity J-fi n 
as a function of the deviation AT of the protocol du- 
ration from the duration T used in Eq. ( 13 1 to calcu- 
late the superadiabatic corrections. These data, shown 
in Fig. ([9]) as black circles, are compared to the results of 
a protocol without the superadiabatic corrections for the 
coupling strength (u/ = uj), while T'(t) follows Eq. (13), 
including the delta-like discontinuities of at the beginning 
and at the end of the protocol (red open squares data of 




FIG. 9: (Color online) Measured fidelity .7-fin of the superadi- 
abatic tangent protocol as a function of the relative deviation 
of the protocol duration from the optimum value T used in 
Eq. ( |13[ ) to calculate oj'(r) and T'(t) (closed black circles), 
and the measured fidelity of a protocol without the superadi- 
abatic corrections for the coupling strength, i.e. u)' — uj (open 
squares). The inset reports a detail of the graph on a linear 
scale. The target T = 1, denoted by the horizontal dashed 
line, is reached at shorter times in the optimized protocol. 
The continuous vertical line shows the minimum time neces- 
sary to achieve T = 1 according to the quantum speed limit. 
Data for uj = 0.5 and T=5.9. 



Fig. (|9j)). This protocol, apart from the discontinuities in 
r'(r), is essentially a tangent protocol. The results are 
summarized in Fig. ^ , which shows clearly that both the 
tangent and the superadiabatic tangent protocols are ex- 
tremely robust with respect to variations of the protocol 
duration. Our study demonstrates that only for large re- 
ductions of the protocol duration does the fidelity drop 
sharply, as otherwise the quantum speed limit would be 
violated. Nevertheless, the tangent protocol reaches a fi- 
nal fidelity J 7 ^ ~ 0.99 for a protocol duration a factor 
two longer than the corresponding superadiabatic proto- 
col. Moreover, the superadiabatic tangent protocol, en- 
suring a perfect following of the adiabatic ground state 
for AT = 0, yields J 7 ^ > 0.99 for an increase in T up to 
100% of its optimum value. 



ite pulse protocol realized such minimum time, T = tt/uj 
given by the quantum speed limit. Finally, a very impor- 
tant indicator is the robustness of the final fidelity respect 
to variations of the sweep parameters. With respect to 
this problem, the superadiabatic protocols, based on a 
completely different approach, have been shown to be 
extremely robust. In fact, while the final fidelity for the 
composite pulse protocol and the RC protocol (the two 
protocols examined in this paper, apart from the supera- 
diabatic ones, requiring the knowledge of the system pa- 
rameters in order to compute the sweeps of ui and T) 
have a non-trivial behavior as a function of the system 
parameters, the superadiabatic protocols fidelity remains 
■^"fin ^ 0.99 for a wide range of parameters. 

In this section the performance indicator is the mini- 
mum time required to reach a final state with J-fi n = 1, 
given some finite experimental resources. 



A. Average or peak coupling strength 

An important question in the use of quantum control 
protocols are the resource requirements. The superadia- 
batic protocols we have investigated require a well chosen 
time dependence T(t) corresponding to the laser detun- 
ing, and in addition they require an increase of the cou- 
pling strength, from uj to uj'(t), which corresponds to an 
increase of the applied laser intensity. On the basis of 
the quantum-speed limit T n w vr/o;, that increase corre- 
sponds to a faster transfer to the target state even with a 
simple Rabi-pulse protocol. Therefore it is important to 
question whether the extra resources should be devoted 
to the application of the superadiabatic protocol, or to 
increasing the performance of a standard protocol. 

The request on the resource increase should be dis- 
cussed in terms of either the peak coupling strength uj pea k 
required for the protocol realization or of the required 
coupling strength averaged over the applied time. We 
introduce the average coupling strength (uj) defined as 



Uj'( T )dT. 



(18) 



Ref. [23] has pointed out that in the case of a time de- 
pendent Hamiltonian the quantum speed limit of Eq. ^ 
should depend on (uj) instead of uj. 



VII. BEST USE OF RESOURCES 

The protocols discussed so far, from generalized LZ 
sweeps to the superadiabatic protocols, can be compared 
in terms of various performance indicators. The first one 
is the capability to reach, at least in theory, a fidelity 
equal to 1 in a finite sweep duration. All the protocols 
examined here, except the standard LZ sweep, have this 
capability. A second indicator is the minimum time re- 
quired to reach Ta n = 1 for a given value of the coupling 
strength uj, constant during the protocol. The compos- 



B. Superadiabatic linear protocol 



The expression for uj'(t) of Eqs. ( 12 ) can be inverted in 



order to derive the initial uj and its peak value during the 
time evolution required to reach the target state within 
a given duration of the protocol T. It turns out that for 
each given duration there exists a unique choice of those 
two parameters leading to a maximum evolution speed. 
This analysis, applied to a range of values, is reported in 
Fig. 10 The dotted line is the prediction of the duration 
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FIG. 10: The red dotted line on top shows the time T required 
for the LZ protocol to reach a 98% final fidelity vs Q. = uj, 
which is time-independent. For the superadiabatic linear pro- 
tocol, the blue continuous line (second from the top) shows the 
time T required to reach 100% final fidelity vs Q = uj pea k, the 
minimized peak value of the time dependent ui'(t), while for 
the blue dashed line (third from top) fi = (oj) . For the supera- 
diabatic tangent protocol, the black dashed-dotted line shows 
the time T required to reach 100% final fidelity vs Q, = uj' , the 
minimized, time-independent coupling constant. The contin- 
uous black line represents the quantum speed limit time T vs 
Q = u>. 



of the LZ protocol required to reach Tfi n = 0.98. The 
other lines report the results for the superadiabatic pro- 
tocols, which have a final fidelity equal to 1. For a given 
duration T of the protocol, we determine the smallest 
values of initial uj and (uj) producing a superadiabatic 
linear transfer within the time T . The continuous blue 
line reporting the dependence of T on the peak value of 
uj' lies below that of the standard LZ protocol. At a given 
peak value, the superadiabatic linear protocol produces a 
perfect target transfer with a 35% decrease in the trans- 
fer time with respect to the LZ protocol with JF^ n = 0.98. 
The dependence of T on (uj), given by the blue dashed 
line, exhibits a change in slope in the intermediate ex- 
plored region, because at small values of uj the linear 
superadiabatic protocol imposes a uj'(t) which differs 
from uj essentially only in the region around r = 0.5 
((uj) 7^ ijjpeak)- The transfer time T for the supera- 
diabatic linear protocol at a given value of (uj) is 50% 
shorter than the time required by the LZ protocol with 
-^fin = 0.98 for (uj) < 1. For (uj) ^> 1 the correction of the 
coupling strength according to the superadiabatic linear 
protocol vanishes (uj' ss uj and (uj) s» uj pea k)- 
In the explored range of the coupling strength, the mini- 
mum time T required by the superadiabatic linear proto- 
col is larger than the quantum speed limit time (continu- 



ous black line in Fig. 10) by a factor between 2 ((uj) sa 5) 
and 10 ((uj) « 0.1). 



C. Superadiabatic tangent protocol 

For the superadiabatic tangent protocol we also ap- 
plied the optimization approach minimizing uj and uj 1 at 
a given T value. Note that for this protocol the coupling 
strength is constant, i.e. uj' = (uj). The results for T vs 
uj' are reported in Fig. [To] (black dashed-dotted line). In- 
terestingly, the optimized protocol reaches the quantum 
speed limit for uj' < 0.1 and deviates at larger uj' . This 
result is quite interesting, since it means that there exists 
a superadiabatic protocol that ensures a perfect follow- 
ing of the instantaneous ground state and approaches the 
quantum speed limit. 



VIII. CONCLUSIONS 

We have explored several quantum protocols transfer- 
ring the initial state of a two-level system into a final or- 
thogonal one. The flexibility and quantum control associ- 
ated with our two level system based on an optical lattice 
BEC configuration has allowed us to measure the final fi- 
delity, the minimum time required to reach the target 
state, and also to determine, in different bases, the tem- 
poral evolution of the wavefunction along the quantum 
trajectory. In addition, we have explored the robustness 
of the final fidelity against a variation of the Hamiltonian 
parameters. We have compared our experimental results 
to numerical simulations based on the Schrodinger evo- 
lution of a single particle. 

Most of our protocols operate with constant uj, hence 
the temporal dependence of F'(t) determines the over- 
all protocol. Figures [2] and [5j where those dependen- 
cies are depicted, shows that the protocol efficiency is 
optimal when T'(t) is nearly constant near the energy 
level anticrossing. Because the adiabatic theorem links 
the tunneling probability at the anticrossing to the time 
derivative of T'(r), a nearly constant value implies that 
the probability of leaving the adiabatic quantum state 
is very small, and the protocol operates with high effi- 
ciency. 

As explored in our previous work |27j . the quantum 
speed limit can be approached using the superadiabatic 
protocol operating in the optimized regime. An impor- 
tant result of the present work is that the quantum speed 
limit is also reached by the power law protocol with large 
a. Therefore a simple modification to the LZ scheme may 
lead to a large improvement in operation speed. We have 
also discussed the important question of the resources to 
be implemented for different protocols in order to reach 
the final state in a given time. 

While the BEC two-level configuration we used in the 
present investigation is not useful as a qubit for quan- 
tum computation, our results can be implemented on 



10 



any qubit configuration, and in particular they could be 
exploited by the trap ion community where very high 
fidelity levels are commonly reached. 
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Appendix A: Two-level BEC Hamiltonian 

The generalized two-level LZ theory is applied to 
study the temporal evolution of ultracold atoms loaded 
into a spatially periodic potential ^ cos (2nx/dL + 4>{t)). 
In the reference frame of the standing wave [35], this 
potential becomes the sum of a periodic potential 
^ cos (2-Kx/di,) and a slowly varying force F in the 
presence of negligible atom-atom interactions, as is in 
the case of our experimental conditions [29] . The char- 
acteristic energy scale of the system is the recoil en- 
ergy defined as E ICC = frw rcc = Tr 2 h 2 /2Mdl. For 87 Rb 
cj roc = 2ir x 3.125 kHz defines the natural units of energy 
huj ICC and time l/u> rcc for our system. 

The lowest energy levels of such a system are given by 
the quadratic dispersion relations of free particles. Mo- 
menta differing by 2p rec (where p rec = hk^) are coupled 
through a ui coupling constant given by 

ho=^, (Al) 

leading to energy bands of the quasimomentum q. The 
dynamics of ultracold atoms in a tilted optical lattice can 
be described by the well-known Wannier-Stark Hamilto- 
nian [3D]. Under the action of the force F, the quasi- 
momentum of a condensate initially prepared at q = 
in the n = energy band scans the lower band in 
an oscillating motion with the Bloch period given by 
TbiocH = 2ith/{FdL)- At the edge of the Brillouin zone, 
around q = n/(2di J ), where an anticrossing of order Vo is 
present, tunneling of the condensate to the n = 1 energy 
band may occur. The full dynamics of the Wannier-Stark 
system can be locally approximated by a simple two-state 
model, with an effective Hamiltonian brought into the 
form of Eq. (fTl) by subtracting the quadratic term in the 
momentum [30]. 

The time dependence of T is induced through a vari- 
ation of the quasimomentum q of the condensate by the 
force F. The time dependence q(t) leads to 

r < f >= 2r »(f -£)• < A2 > 



where ±ro represent the initial and final values of F. 

There are some limiting cases and experimental pa- 
rameters for which the simplified two-state model is not 
a good approximation for our system. The discrepancy 
is considerable for lattice depths very large compared to 
the energy scale E Tec of the system, for which the gap be- 
tween energy bands increases leading to quasi-flat bands 
and localized eigenstates. Therefore, several momentum 
eigenstates contribute with a non-negligible amount to 
the lowest energy eigenstate, and one would need to take 
into account more components in the Hamiltonian ma- 
trix. 

While the final fidelity values are obtained performing 
a projective measurement after accelerating the lattice 
for a full Bloch period, the F(t) fidelity is obtained by 
performing a projective measurement at a time t smaller 
than T Bioch- 

Appendix B: Experimental implementation 

In our experiments we realized an effective two-level 
system using ultracold atoms, forming a Bose-Einstein 
condensate, in an optical lattice [29] subjected to a 
time-dependent force. Initially, we created Bose-Einstein 
condensates of 5 x 10 4 87 Rb atoms inside an optical 
dipole trap (mean trap frequency around 80 Hz). A 
one-dimensional optical lattice potential created by two 
counter-propagating, linearly polarized Gaussian beams 
was then superposed on the Bose-Einstein condensate by 
ramping up the power in the lattice beams to a final Vq 
value in the range 1 to 5 E rcc . The wavelength of the 
lattice beams was Xl = 842 nm, leading to a lattice con- 
stant d L = A/2 = 421 nm. 

The experimental protocols are carried out using the 
techniques previously developed by us and described in 
detail in j^H [3D]. The force is experimentally imple- 
mented by accelerating the optical lattice [33]. A small 
frequency offset Au(t) between the two beams was intro- 
duced through the acousto-optic modulators in the setup, 
which allowed us to accelerate the lattice in a controlled 
fashion and, hence, to control the quasi-momentum q(t) 
of the condensate. The time dependence of to is controlled 
through the power of the lattice beams which determines 
V . 

The final state is detected by measuring the atomic mo- 
mentum distributions in time-of-flight images. Our BEC 
has an initial width in momentum space much smaller 
than the width of the first Brillouin zone. This enables 
us to observe the full dynamics for single or multiple 
crossings, the only limitation being the initial momen- 
tum width of the condensates and nonlinear effects. 

The noise in our imaging system does not allow us 
to reliably to measure fidelities J- > 0.98. In addition, 
the non-adiabaticity of the preparation and measurement 
protocols contributes with an infidelity on the order of 
0.01. 



11 



[1] M.V. Berry, Ann. N.Y. Acad. Sci. 755, 303-317 (1995). 

[2] L. Landau, Phys. Z. Sowjetunion 2, 46 (1932). 

[3] C. Zener, Proc. R. Soc. A 137, 696 (1932). 

[4] The evolution of a two-level system was investigated at 

the same time in [5] , and the tunneling is also identified 

as Landau-Zener-Majorana tunneling. 
[5] E. Majorana, Nuovo Cimento 9, 4350 (1932). 
[6] I. Walmsley, and H. Rabitz, Phys. Today 56, 43 (2003). 
[7] S.A. Rice, and M. Zhao, Optical Control of Molecular Dy- 
namics (Wiley, 2000). 
[8] B. M. Garraway, and K.A. Suominen, Rep. Prog. Phys. 

58 365-419 (1995). 
[9] M.H. Levitt, M.H. Composite Pulses, in: Encyclopedia of 

Nuclear Magnetic Resonance, eds. D.M. Grant, and R.K. 

Harris (Wiley, 1996). 
[10] J. Roland, and N.J. Cerf, Phys. Rev. A 65, 042308 

(2002). 

[11] D. A. Garanin and R. Schilling, Europhys. Lett. 59, 7 
(2002). 

[12] D. A. Garanin and R. Schilling, Phys. Rev. B 66, 174438 
(2002). 

[13] R.Lim, and M.V. Berry, J. Phys. A: Math. Gen. 24, 3255 
(1991). 

[14] M.V. Berry, J. Phys. A: Math. Theor. 42, 365303 (2009). 
[15] M. Demirplak, and S.A. Rice, J. Phys. Chem. A 107, 
9937 (2003) 

[16] M. Demirplak, and S.A. Rice, J. Chem. Phys. 129, 
154111 (2008). 

[17] J.G. Muga, X. Chen, A. Ruschhaupt, and D. Guery- 
Odelin, J. Phys. B: At. Mol. Opt. Phys. 42 241001 (2009). 

[18] A. Ruschhaupt, X. Chen, D. Alonso, and J.G. Muga, 
New J. Phys. 14 093040 (2012). 

[19] Levitin, L.B. Int. J Theor. Phys. 21, 299 (1982). 

[20] Bhattacharyya, K. J. Phys. A:Math. Gen. 16, 2993 
(1983). 

[21] Peres, A. Quantum, Theory: Concepts and Methods 

(Kluwer, 1993). 
[22] V. Giovannetti, S. Lloyd, and L. Maccone, Phys. Rev. A, 

67, 052109 (2003). 
[23] T. Caneva, M. Murphy, T. Calarco, R. Fazio, S. Mon- 

tangero, V. Giovannetti, G. E. Santoro, Phys. Rev. Lett. 

103, 240501 (2009). 
[24] A. Mitra, A. Ghosh, R. Das, A. Patel, and A. Kumar, J. 

Mag. Res. 177 285298 (2005). 
[25] T. Lu, X. Miao, and H. Metcalf, Phys. Rev. A 75, 063422 

(2007). 

[26] A. S. Mellish, G. Duffy, C. McKenzie, R. Geursen, and 
A. C. Wilson, Phys, Rev. A 68, 051601(R) (2003). 

[27] M.G. Bason, M. Viteau, N. Malossi, P. Huillery, E. Ari- 
mondo, D. Ciampini, R. Fazio, V. Giovannetti, R. Man- 
nella and O. Morsch, Nat. Phys. 8, 145 (2012). 

[28] J.-F. Schaff, P. Capuzzi, G. Labeyrie, and P. Vignolo, 
New J. Phys. 13, 113017 (2011). 

[29] A. Zenesini, H. Lignier, G. Tayebirad, J. Radogostowicz, 
D. Ciampini, R. Mannella, S. Wimberger, O. Morsch, 
and E. Arimondo Phys. Rev. Lett. 103, 090403 (2009). 

[30] G. Tayebirad, A. Zenesini, D. Ciampini, R. Mannella, 
O. Morsch, E. Arimondo, N. Lorch, and S. Wimberger, 
Phys. Rev. A 82, 013633 (2010). 

[31] M.V. Berry, Proc. Roy. Soc. Lond. A 429, 61 (1990). 

[32] K.W. Madison, C.F. Bharucha, P.R. Morrow, 



S.R.Wilkinson, Q. Niu, B. Sundaram, M.G. Raizen, 
Appl. Phys. B 65, 693700 (1997). 
[33] B. P. Anderson and M. A. Kasevich, Science 282, 1686 
(1998). 



